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ATTRACTIVITY, DEGENERACY AND CODIMENSION OF 
A TYPICAL SINGULARITY IN 3D PIECEWISE SMOOTH 

VECTOR FIELDS 

TIAGO CARVALHO^ AND MARCO A. TEIXEIRA^ 


Abstract. We address the problem of understanding the dynamics 
around typical singular points of 3D piecewise smooth vector fields. A 
model Zq in 3D presenting a T-singularity is considered and a complete 
picture of its dynamics is obtained in the following way: (i) Zq has an 
invariant plane ttq filled up with periodic orbits (this means that the 
restriction ZoLo i® ^ center around the singularity), (it) All trajectories 
of Zo converge to the surface ttq, and such attraction occurs in a very 
non-usual and amazing way, (iii) given an arbitrary integer k > 0 then 
Zq can be approximated by vro-invariant piecewise smooth vector fields 
Ze such that the restriction Ze\-nQ has exactly fc-hyperbolic limit cycles, 
(iv) the origin can be chosen as an asymptotic stable equilibrium of Ze 
when fc = 0, and finally, (v) Zq has infinite codimension in the set of all 
3D piecewise smooth vector fields. 


1. Introduction 

Our interest in this paper is to study some qualitative aspects of piecewise 
smooth vector fields (PSVFs for short) in 3D. Roughly speaking, as stated 
in [5], a PSVF Z = {X,Y) is a pair of C""-vector fields X and Y (both 
defined on M^), in such a way that just their restrictions to some regions 
(half-spaces) separated by a codimension one surface S (called switching 
manifold) are considered. In Section [2] we give a precise dehnition. 

In this context, some of the points with richer dynamics are those ones 
where the trajectories of X and/or Y are tangent to S. These points are 
called tangential singularities. The most known tangential singularity of a 
smooth system X is the fold singularity (also called fold point), which is 
characterized by the quadratic contact of an orbit of X with S. A fold 
point p can be visible or invisible. It is visible for X if the A-trajectory 
passing though p remains in the same side where X is dehned, otherwise it 
is invisible. In 3D, generically there exists a curve of tangential singularities 
Sx C S passing through p (the same for Y). 

If p G S is a fold singularity of both systems X and Y then it is called a 
two-fold singularity. This singularity is a prototypical model in the generic 
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classification of singularities in PSVFs. As pointed out in [7], a two-fold 
singularity is an important organizing centre because it brings together all 
of the basic forms of dynamics possible in a PSVF. There are many dis¬ 
tinct topological types of two-fold singularities and the most interesting of 
them is the so called T-singularity. We say that p is a T-singularity (or 
Teixeira-singularity — due to the pioneering work m — or invisible two¬ 
fold singularity) for Z = (A, Y) if p is an invisible fold point of both X and 
Y and Sx meets Sy transversally at p (see Figured]). The interested reader 
can see more details about the T-singularity in [a El El El El da da eq). 



Figure 1. T-singularity. 

Our primary concern in this article is to show the existence of a infinite 
codimension two-fold singularity. In order to do it, we study smooth non¬ 
linear perturbations of a specific model of 3D degenerate T-singularity and 
a complete picture of its amazing dynamics is exhibited. It is worth noting 
that some interesting bifurcations of such T-singularity are exhibited. 

A formal codimension study of a singularity must contemplate the mini¬ 
mal number of parameters necessary in order to obtain all topological types 
of dynamical systems around the dynamical system presenting the singular¬ 
ity. In [5] the authors perform this study for a planar two-fold singularity 
of PSVFs. 

Some ideas and constructions presented in [5] and [6] are adapted in our 
approach. As consequence, we are able to produce an arbitrary number of 
topological types of PSVFs in a small neighborhood of a PSVF presenting 
a T-singularity. So, this T-singularity has infinite codimension. 

Regardless of, this work fits into a general program for understanding the 
dynamics of higher dimension vector fields expressed by: 

(1) ii = F{u) + sgn{f{u))K 

where u = {ui,U 2 , ■ ■ ■ ,Un) G M”, K is a constant re-dimensional vector, 
/ : M” ^ M is a smooth mapping and F : M” —?■ M"", in general, is not 
smooth in the switching manifold S = {/~^(0)}. 

Systems like ([I]) are wide used in applications and appears in models of 
Mechanics Engineering (see dEolEHl), Electric Engineering (see [a ED), 
Biologic/Control Theory/Economics (with sudden external influences, see 
among others. In fact, every system susceptible to on-off operations 
are modelled by systems like ([T|) which imposes an interdisciplinary aspect 
on this theory. 
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The sgn function produces, in general, lost of differentiability of ([I]) when 
the trajectories pass throught S. So, m is a non-smooth (or piecewise 
smooth) system. As also happens for smooth systems, some choices on 
the functions at the right-side of ([T]) generate non stable behavior. There 
are in the current literature a huge variety of papers dealing with stability 
conditions of models like 0 (see [ElEo] among others). 

It is worth mentioning that in [T] Anosov proves the asymptotic stability 
of the origin of ([H) when T is a linear vector field. After, Kiipper-Kunze- 
Hosham found invariant varieties of ([T]) when K = 0. In fact, they show the 
existence of invariant varieties (cones) and prove that the trajectories of (fTj) 
twist (curl up) until it (see [T^ [T6l [T7] j. 

In our approach the following 3D system is considered: 


Zo{x, y, z) = (x, y, z) = F{x, y, z) -|- sgn{z)K, 
where f{x, y, z) = 2 ;, 


F{x,y,z) 


1 

2 


-1 -I 0 
-1 -1 0 
1 -1 0 


\+-sgn{z)\ 


-1 

-1 

-1 


-1 0 

-1 0 

-1 0 


and 


or equivalently: 



X 

y 

z 


( 2 ) 

Zo{x,y,z) 


X{x, y, z) = (-1 - {x + y),l - {x + y), -y) if 2 > 0, 
Yo{x,y,z) = {l,-l,x) if 2 < 0. 


In fact, as exhibited in [Ml IE], we also identify an invariant manifold in 
our model. Morevoer, the trajectories converge to it in a very unusual way. 
The richness of our model is revealed after a suitable perturbation of it. In 
such case we conserve the invariant manifold and identify the birth of an 
arbitrary number of limit cycles. Also, we obtain the asymptotic stability 
of the origin (as in [I])- 

Now we state the main results of the paper. 


Theorem A. Let Zq be given by ([2]) and ttq the plane {y + x = 0}. For each 
integer k > 0, there exists a one-parameter family of TTQ-invariant PSVFs 
Z^ satisfying: 

(a) Zg — >• Zq when e —>■ 0; 

(b) Ze has exactly k hyperbolic limit cycles in a neighborhood of the 
origin. The same holds for k = 00 and, 

(c) All trajectories of Zq and Z^ converge to ttq. 

(d) When e > 0 and k = 0, the origin is an asymptotic stable equilibrium 
point for Zg. 
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An immediate consequence of Theorem A is; 

Theorem B. The PSVF Zq, given by has infinite codimension. 

The paper is organized as follows. In Section [2] we introduce the termi¬ 
nology, some definitions and the basic theory about PSVFs. In Section[3]we 
present properties towards the understanding of the phase portrait of ([2]). 
In Section 0] suitable perturbations of ([2|) are considered and the birth of 
limit cycles are explicitly exhibited. In Section [5] we prove the main results. 
In Section [6] we made a brief conclusion about the results in this paper and 
in Section [7] we picture some numerical analysis and the phase portrait of 
([2j) around the origin. 


2. Preliminaries 


In this section we give a brief review of the theory. 

Let V be an arbitrarily small neighborhood of 0 G M^. We consider a 
codimension one manifold S of given by S = /“^(O), where / : y —)• ffi 
is a smooth function having 0 G M as a regular value (i.e. V/(p) 0, for 

any p G /“^(O)). We call S the switching manifold that is the separating 
boundary of the regions S"'' = {q ^V \ f{q) > 0} and = {q \ f{q) < 
0}. Throughout paper we assume that S = /“^(O), where f{x,y,z) = z. 

Designate by x the space of C^-vector fields onV C endowed with the 
C^-topology, with r > 1 large enough for our purposes. Call TT the space of 
vector fields Z : 1/ —)• such that 


(3) 


Z{x,y,z) 


X{x,y,z), for (x,y, 2 ;)GS+, 
y{x,y,z), for {x,y,z)eT.-, 


where X = {Xi,X 2 ,Xfi),Y = (Yi,y 2 ;P 3 ) £ X- We endow TT with the 
product topology. The trajectories of Z are solutions of g = Z{q) and 
we will accept it to be multi-valued in points of S. The basic results of 
differential equations, in this context, were stated in m- 

On S we distinguish the following regions: 

• Crossing Region: = {p G S| A 3 (p) • Y^^p) > 0}. Moreover, we 

denote = {p G S | X 3 {p) > 0 ,Y3{p) > 0} and = {p G S | Xfip) < 
0,Y3{p) < 0 }. 

• Sliding Region: S® = {p G S | ^ 3 (p) < 0, T 3 (p) > 0}. 

• Escaping Region: S® = {p G S | X 3 {p) > 0, l 3 (p) < 0}. 

When q G S*US®, following the Filippov’s convention, the sliding vector 
field associated to Z G D’' is the vector held Z® tangent to U and 
expressed in coordinates as 

" yiX^fiq), (X 2 E 3 - y2X3){q), 0 ). 

Associated to (|4j) there exists the planar normalized sliding vector field 


(5) Z^{q) = {{X 1 Y 3 - yiX^fiq), (^ 2^3 - >2X3)(g)). 
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Remark 1. Note that, if q € then X 3 {q) < 0 and Y^{q) > 0. So, 
(I 3 — X'i){q) > 0 and therefore, and are topologically equivalent in S® 
since they have the same orientation and can be (J-extended to the closure 
S* ofTi^. Ifq&T,^ then Z^ and Z^ have opposite orientation. 

In this context, a rich dynamics occurs on those points p £ T, such that 
Xsi^p) - ysip) = 0, called tangential singularities of Z (i.e., the trajectory 
through p is tangent to S). 

For practical purposes, the contact between the smooth vector field X 
and the switching manifold S = /“^(O) is characterized by the expression 
X.f{p) = {Vf{p),X{p)) = 0 where (.,.) is the usual inner product in In 
this way, we say that a point p G S is a fold point of X if X.f{p) = 0 but 
X'^.fip) / 0, where X^.f{p) = (XX^~^. f (p), X{p)'j for i > 2. Moreover, 
p € S is a visible (respectively invisible) fold point of X if X.f{p) = 0 
and X‘^.f{p) > 0 (respectively X‘^.f{p) < 0). In addition, a tangential 
singularity q is singular if g is a invisible tangency for both X and Y. On 
the other hand, a tangential singularity q is regular if it is not singular. 

Call Sx (resp. Sy) the set of all tangential singularities of X (resp. 
y). In ?)D, a point of S at which two curves of fold singularities Sx and 
Sy meet is called a two-fold singularity. This singularity is a prototypical 
model in the generic classification of singularities in PSVFs. As pointed 
out in [7], a two-fold singularity is an important organizing centre because 
it brings together all of the basic forms of dynamics possible in a PSVF. 
There are many distinct topological types of two-fold singularities and the 
most interesting of them is the so called T-singularity. We say that p is 
a T-singularity (or Teixeira-singularity or invisible two-fold singularity) for 
Z = [X, Y) if p is an invisible fold point of both X and Y and the intersection 
of Sx and Sy is transversal at p (see Figure HI)- It is easy to check that 
in the model ([2]) the origin is a T-singularity. In this article, we study 
smooth nonlinear perturbations of this model and a complete picture of its 
dynamics is exhibited. The interested reader can see more details about the 
T-singularity in O [3l EJ El El E2l EHl |20] . 

It is worth to say that some constructions and ideas of [5] and [ 6 ] are very 
useful in our approach. 

Now we establish a convention on the trajectories of orbit-solutions of a 
PSVF. 

Definition 1. The local trajectory (orbit) (pzi^iP) of a PSVF given by 
through p £V is defined as follows: 

• For p G S+\S = {(? G V I z > 0} and p G S“\S = {(? G V | z < 0} 
the trajectory is given by (l)z{t,p) = 4>x{t,p) and (fzit^p) = 4>Y{t,p) 
respectively. 

• For p G and taking the origin of time at p, the trajectory is 
defined as (f)z{t,p) = 4>y{t,p) for t < 0 and (j)z{t,p) = 4>x{t,p) for 
t > 0. For the case p G the definition is the same reversing 
time. 
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• Forp G S® and taking the origin of time atp, the trajectory is defined 
as 4>z{t,p) = (f)z'^{t,p) for t < 0 and 4>z{t,p) is either (f>x{t,p) or 
fivit-iP) or 4>2^{t,p) for t > 0. For p G S® the definition is the same 
reversing time. 

• For p a regular tangency point and taking the origin of time at p, the 
trajectory is defined as fizit^p) = 4>i{t,p) for t < 0 and (j)z{t,p) = 
4 > 2 {t,p) for t > 0, where each (pi, (p 2 is either cpx or cpy or (pz^- 

• For p a singular tangency point (pz{t,p) = p for all i G K. 

Definition 2. The orbit (trajectory) of a point p £ V is the set ^{p) = 
{(pzit,p) : t G M} obtained by the concatenation of local trajectories. 

Consider 0 / p G It is easy to see that there exists a time ti{p) > 0, 
called time, such that the forward trajectory of X passing through p 

at t = 0 return to S after ti (p) . We define the half return map associated to 
X by px{p) = 4>x{fi{p),p) = Pi G S. When pi G let t 2 (pi) > 0 be the 
Wfly time of the trajectory of Y passing through pi. Define the half return 
map associated to Y by Py{pi) = f>Y {t 2 {pi) ,Pi) £ S. The involution px 
(resp. py) is such that Fixfpx) = Sx (resp. Fix{pY) = Sy). The first 
return map associated to Z = {X,Y) is defined by the composition of these 
involutions, i.e., 

(6) pz{p) = py o px{p) = 0 y(i 2 (pi), (px{h{p),p)) 

or the reverse, applying first the flow of Y and after the flow of X. See 
Figure [2] and details in m- 



Figure 2. First Return Map. 

The mapping pz is an important object in order to study the behavior 
of Z around a T-singularity. The proofs of the main results require a de¬ 
tailed analysis of the return map and different domains (departure regions), 
according to the Filippov’s decomposition of S, must be considered. For 
example, if g G then the iteration of both mappings py o px{q) and 
Px ° Ty{q) must be considered. 

3. Properties of System Zq given by ([2}) 

In this section we describe some important features about the PSVF given 
by ([2]). In Subsection 13.11 we describe the sliding vector field associated to 
([2|). In Subsection 13.21 we analyze the first return map associated to it. In 
Subsection 13.31 we analyze the way in which the trajectories converge to a 
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limit set. This last analysis permits us to detect the behavior of certain 
invariant sets. 

3.1. The sliding vector field associated to ([2]). Using Equation (j5|). 
associated to ([2|) we have the normalized sliding vector field 

(7) Z^(x, y, z) = ((-1 - (x + y))x + y, (1 - (x + y))x - y, 0). 

Let us understand the phase portrait of Q. 

Proposition 3. The normalized sliding veetor field Zq, given by 0, has a 
saddle-node at the origin. 

Proof. Identify S with the xy-plane. Consider the change of variables (rt, u) = 
(x + y, X — y). Then ([7]) can be re-written in the form {u,v) = {—{u + 
v)u, —2v). This last system has the origin as a unique equilibrium with 
eigenvectors vi = (1,0), V 2 = (0,1) associated to the eigenvalues Ai = 0, 
A 2 = —2 respectively. The phase portrait is pictured at Figure [3l 



Figure 3. Phase Portrait of the normalized sliding vector field Zq. 


□ 

Remark 2. Since d?]) has a saddle-node at the origin, by Remark\^we must 
reverse the orientation in So, we conclude that the sliding vector field 
associated to Zq has the phase portrait shown at Figure Note that the 
straight line y -\- x = 0 in Ti, where Yq — Xq = 0 in ([4|), is composed only 
by equilibrium points of the sliding vector field associated to Zq. Moreover, 
exeept by the stable invariant manifold, the trajeetories of the sliding veetor 
field associated to Zq depart from S®. 

3.2. The first return map associated to ([2]). In order to exhibit the first 
return map associated to Zq = {X,Yq), given by ([2]), we have to write the 
expressions of the half return maps tpx and ipYo ■ A straightforward calcula¬ 
tion shows that the trajectories 4>x and fivo are parameterized, respectively, 
by 

( 8 ) (pxit) = {-t-\-kie~‘^^-\-k2,t-\-kie~‘^* -k2, /2 + kie~'^^/2 + k2t + kfi), 
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Figure 4. Phase Portrait of the sliding vector field associated to Zg. 


and 

(9) 4>Yo(t) = {t + li, —t + ^2) i^/2 + ht + lf/2 + /s), 

Consider the planes = {(x, y,z) G V \ y = —x + A:}, with A: e M. 

Proposition 4. Given an arbitrary point (xo,yO)0) G then 

2/0,0) = (-A2 + ((a:o+yo)/2)e"2*2 _^(xo-2/o)/2, t2 + {{xQ+yQ) /2)e-‘^^^ - 
(xo — 2/o)/2,0), where the X-fly time t 2 > 0 is given implicitly by 

(10) - A2/2 + ((xo + yo)/ 4 )e “^‘2 _ yQ)/2)t2 - (xq + yo)/4 = 0. 

In particular, (j)x{T^QGTj'^) C (tto n S+) and ipx{xQ,—XQ,i)) = (— xo,xo,0). 

Proof. Considering the initial condition (xo,2/o,0) and ([8|), let t 2 > 0 be the 
first time such that 

-t^/2 + ((xo + yo)/4)e“^‘ + ((xq - yo)/2)t - (xq + yo)/4 = 0. 

Using ([8]) it is easy to see that 

Px{xo,yo,0) = 

{-t 2 + {{xo+yo)/‘2)e~^^^ + {xo-yo)/‘^,t 2 + i{xo+yo)/‘^)e~^^^ -{xo-yo)/2,0). 
In particular, for pQ = (xo,2/o,0) = (xq,— xo,0) G tto n S, we obtain that 
the first two coordinates of 4>x{t,Po) are x(t) = —A + xq and y{t) = t — 
Xq = —x{t). Moreover, in this case we can solve explicitly the equation 
-A^/2 + ((xo + yo)/2)e“2*V2 + ((xq - yo)/2)t2 - (xq + yo)/4 = 0 and we 
obtain t 2 = 2xo. So, <y9j\:(xo, —xq, 0) = (—xo,xo,0). This concludes the 
proof. □ 

Proposition 5. Given an arbitrary point (xo,yo,0) G then 

PYo{xo,yo,0) = (-xo,yo + 2xo,0). In particular, (/>yo(7rfc n S") C (TTfcnS"). 

Proof. Considering the initial condition (xq, yo, 0) and ([9]), in order to deter¬ 
mine (fYo it is enough to obtain the first time Ai > 0 such that A^/2-|-xoAi = 0. 
So Ai = —2xo and (/?yg(xo, yo, 0) = (—xq, yo + 2xo, 0). In particular, for 
Po = {xo,yo,0) = (xo,—Xo -I- A:,0) G vr^ n the first two coordinates of 
(pY{t,Po) are x(t) = t + XQ and y{t) = —t — XQ + k = —{t+xo) + k = —x(t) + k. 
This concludes the proof. □ 
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Proposition 6. The plane vro = {(x, y, z) G P | x + y = 0} is Zo-invariant. 
Moreover, ZqItto is a center. 

Proof. By Propositions [4] and [5] we get that ttq is invariant by the flow of 
Zq. In order to see that Zq has a center at tto it is enough to see that 

TZoixo,-xo,0) = ipYo o ipxixo,-xo,0) = ipYoi-xo,xo,0) = (xo,-xo,0). 

□ 


Now we will prove that tto is a hyperbolic global attractor for the trajec¬ 
tories of Zq. 

Proposition 7. Let tq be the straight line given by vq = vro H S. Given 
{xq, yo,0) G then 

d{TZo{xo,yo,0),ro) < d((xo, yo, 0), ro). 

This means that the trajectories of Zq are converging to tq. Moreover, 

Tzo{xo,yo,t)) = {-Xn,Xn +{xo + yo)e~^d2 +-+h 

where 4*^ is the fly time necessary to the X-trajectory by ip'^^^{xo,yo,0) re¬ 
turns to S and Xn = —4”^ + ((xo + yo)/ 2 )e “^(*2 ^+...+4 4 _ Xn-i — (xq -|- 

yo)(,-2(4^)+...+4"-^4)/2. 


Proof. By Proposition [4] we obtain 

Tx{xo,yo,0) = 

{-t 2 + i{xo+yo)/2)e~‘^^^ + {xo-yo)/ 2 ,t 2 + {ixo+yo)/2)e~‘^^^-{xo-yo)/2,0), 
where t 2 > 0 is given implicitly by 

-4/2 -t ((xo -h yo)/2)e"^*V2 + {{xq - yo)/ 2 )t 2 - {xo + yo)/4 = 0. 

By Proposition m 
( 11 ) 

TZoixo,yo,0) = {ipY o ipx)ixo,yo,0) = 

= (^t2 - ((£0+M)e-2i2 _ (£0^^ 0Eo+M)e-2i2 + (£2^, q) . 

Then we get, 

\/2 y/2 

d{pzo{xo,yo,0),ro) = —(xo + yo)e“^*2 < — (xq + yo) = d((xo,yo,0),ro). 
In order to obtain that 

V 2 zn(^o,yo, 0 ) = (-x„,x„-h (xo+ yo)e“ ^ 2 +-+2 i^o), 


with 

^ _ (n) (xQ + yo) 2 (tf^+...+tir^) ^ (xo + yo)(e-2(4^^+-+4” 

■^n — ^2 ' 2 ^ •^n—1 2 ’ 

it is enough to use n times Propositions 0] and O □ 
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Proposition 8. Let po = (xo,yO)0) G tiq. The first return map pz{po) 
has eigenvalues pi = 1 and /i2 = 1 — 4 xo + 8 xq + 0 {xq) and eigenvectors 
ui = (-1,1) and U2 = (-(1 - 2xo + 2xo)/(l - 4 xo + 6x9) + 0 {xl),l), 
respectively. 

Proof. Since we can not solve Equation (jlOl) we consider the second order 
series expansion of So we are able to solve (1101) and obtain the X-?^y 

time of pq. Moreover, a straightforward calculation shows that the diffeo- 
morphism pz{po) has eigenvalues pi = \ and /i 2 = 1 — 4xo + 8 xq and 
eigenvectors ui = (— 1 , 1 ) and U 2 = (—(1 — 2xo + 2xo)/(l — 4xo + 6 xq ), 1 ), 
respectively. □ 

Proposition 9. All trajectories of Zq, given by ([2]), converge to the plane 
TTo. 

Proof. First note that each q € V hits T, for some positive time. Using 
Remark [2l the position of the eigenspaces of dH) and Propositions [7] and [8] 
it is easy to see that given p G S = S® U U the trajectories of Zq by p 
converge to tto. □ 

3.3. The convergence of the trajectories. Now we picture the scenario 
describing the asymptotic behavior of Zq in V. As we said above, the plane 
vTo is an attractor for fizit^p) and separates V in two open regions = 
{(x, y,z) G U I X + y > 0} and V~ = {(x, y,z) G U | x + y < 0}. 

Since all points of V hit S, in order to determine the asymptotic behavior 
of a trajectory it is enough to take this trajectory departing just from points 
in S. Consider the partition of S = U U U U U U rjj" U 
r(7 U 5^ U 53 ^ U U 5^ U 0}, where n U+, n U", 

n V-, n U+, r+ = ro n r^ =ron Sy = 

SvnV-,Sf^ = Sxnv-, = 5rnu+, 5^ = 5xnu+ and 0 = {(0,0,0)}. 
See Figure 0 



Consider the following steps: 

(1) First let us analyze the trajectories in tq. 

(l.i) If p = 0 then, according to the fifth bullet of Definition [U 
fiz{t,p) = 0 for all t G M. 
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(1.11) If p G Tq U then, according to Proposition [6l (j)z{p) 
became restrict to vro and describes a periodic orbit around the ori¬ 
gin. 

(2) Let us analyze the trajectories in H S. 

(2.1) If p G then, according to Proposition [3l 4'zit,p) —)• 0 
when t —)• -boo. By the fourth bullet of Dehnition[Tl the same holds 
when p G Sy U . 

(2.11) If p = {xp,yp,0) G then we can put yp = —Mxp, 
where 0 < M < 1. With this, we are able to explicitly calcu¬ 
late the successive returns p^{p) = (n”(p), II^(p), 0), with n > 1. 
A straightforward calculation shows that 112(0) is positive when 
2n/{2n + 1) > M. This implies that after the hrst integer N such 
that 2N/{2N + 1) > M we obtain n 2 (p) > 0 and so W^{p) belongs 
to S®, where y > 0. By the analysis done in (2.i), (j)z{t,p) —)• 0 
when t —>■ -boo. This means that the competition between the two 
attractors ttq and is won by 

(2.111) If p = {xp,yp,0) G then cpzip,—“^Xp) G E^+ and so, 
repeat the analysis of item (2.ii). 

(3) Let us analyze the trajectories in P n E. 

(3.1) If p G Eg"'' U Sy then we can put yp = —Rxp, where 0 < 
R < 1. With this, we are able to explicitly calculate the successive 
returns <p^{p) = (n”(p), II^ (p), 0), with n > 1. A straightforward 
calculation shows that p^(p) ^ oo. So, (pz{t,p) —>• oo when t —)• 
-boo. 

(3.11) If p G Eg U then (j)z{p, —2xp) G Eg"'' and so, repeat 
the analysis of the previous bullet. 

(3.111) If p G E*^ then, according to the third bullet of Definition 

[U there are three choices for (j)z{t,p). When (j)z{t,p) = 4>x(t,p) or 
when (pzip, t) = <pY{t,p) the trajectory hits Eg“ or Eg+j respectively. 
In both cases we use the previous two items. When (f>z{t) = {t,p) 

then, by Remark [21 there exists an invariant manifold in E® converg¬ 
ing to the origin and the other points in E® converge to SyUS^ (and 
so, we use the previous two items). This means that the competition 
between the attractor ttq and the repulsiveness of Z in p G Eg"*" is 
won by the second one. 

Remark 3. As consequence of the previous analysis we are able to say that 
any Z-trajectory, with Z given by ([2]), through p G V either describes a 
periodic orbit or converges to a stationary point or leaves any ball around 
the origin, according to the initial position of p. 


4. Properties of an oriented perturbation of System Zq 
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4.1. Auxiliary results. In what follows, /i : M —)• M will denote the C°°- 
function given by 


h{w) 


0, if u) < 0; 

if n; > 0. 


Consider the function Fi{x, y) :M? x M, where either p = f oic p = i, such 
that 

(12) F/ (x, y) = —eh{x)h{—y){e — x){2£ — x)... [ke — x) 


with fc € N, and 

(13) F^(x,y) = h{x)h{—y)sm{Tr£‘^/x). 

Lemma 10. Consider the function Fi{x,y) given by (jl2p . 

(i) If £ < 0 then fI does not have roots in (0,+oo) x {y}. 

(ii) Ife > 0 then fI has exactly k roots in (0, +oo) x {y}, these roots 
are {(e, y), (2e, y),... ,{k£,y)}. 

(hi) forj e {1,2, ...,A:}. 

It means that such partial derivative at (je, y) is positive for j odd 
and negative for j even. 

Proof. When x > 0, by a straightforward calculation Fi{x,y) = 0 if, and 
only if, (e — x)(2e — x)... (fee — x) = 0. So, the roots of f/ (x, y) in (0, +oo) 
are e, 2e,... , k£. Moreover, 

^^(x,y) = -eh{-y)-^{{j£-x)H{x)^ = -eh{-y){{j£-x)^{x)-H{x)^ , 

where F(x) = f/{ x,y)/{eh{—y){j£ — x)). So, 
f)TpS 

^^(je,y) = eh{-y)H{j£) = 

= eh{-y)£^h{j£){l - j)... ((j - 1) - j)((j + 1) - j) ■ ■ • (^ - j) 

= -eh{-y)£^h{j£){-iy (^{j - 1)... (j - {j - 1))) (^((j + 1) - j)... {k - j)) 
= -eh{-y){-iy£^h{j£){k - j)l{j - 1)! 


This proves items (ii) and (hi). Item (i) follows immediately. □ 

Lemma 11. Consider the function F^{x,y) given by (1131) . For £ y 0 
the function Ft has infinitely many roots in (0,e^) x |yj, these roots are 
{(e^ y), (eV2, y), (eV3, y), • • • } and 

f) 

^(e^/j,y) = h{-y){-iy {-7rf /£'^)h{£'^ /j) forj e {1,2,3,...}. 

It means that such derivative at {£^lj,y) is positive for j odd and negative 
for j even. 
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Proof. When x > 0, by a straightforward calculation F^(x,y) = 0 if, and 
only if, sin(7re^/x) = 0. So, the roots of F^{x,y) in (0,e^) x {y} are 
(e^ y), (eV2, y), (eV3, y), • • • • Moreover, 

f) 

-;^{x, y) = h{—y)[h'{x) sin(7re^/x) + h{x) cos(7r£^/x)(—7re^/x^)]. 
ox 

So, 

-^(eVy>y) = h{-y)[h'{e^/j)s\n{-Ke^/j)+ 

+ /i(£Vj)cos(7re2/j)(-7rjVe^)] 

= h{-y){-iy{-TTf/e'^)h{e'^/j). 

□ 


Consider Zq given by ([2]) and 


(14) 

Zf{x,y,z) 


X{x, y, z) = (-1 - (x + y), 1 - (x + y), -y) if z > 0, 
Yi{x,y,z) = (l,-l,x + ^(x,y)) if z < 0. 


Remark 4. Take = Zf, with p = i,f. It is easy to see that Zg —)• Zq 
w/ien e ^ 0. 


Associated to ()14j) we have the normalized sliding vector field given by 

Zp,e{x,y,z) = 

/ dF^ dF'^ \ 

(^(-1 - (x + y))(x + -^{x, y)) + y, (1 - (x + y))(x + -^{x, y)) - y, Oj. 

A straightforward calculation shows that the trajectory (fyp of Yf given 
in da are parameterized by 

(15) (^) ~ {t Y li, —t Y l2,t^/2 Y lit Y Ff (t Y li, —t Y I 2 ) + ^ 3 )- 

Proposition 12. Given an arbitrary point (xo,yo,0) G then 

</?yp(xo,yo,0) = {ti Y xo,—ti + yo,0), where the Yf-fly time ti > D is 
given implicitly by t\/2 + xoti + Ff{ti Y xq, —ti + yo) = 0. In particular, 
(pYfi^k n S") C (vTfc n S"). 

Proof. Let ti > 0 the first time such that /2+xoti + P/(ti +xo, —ti +yo) = 
0. Using (fTCI) it is easy to see that (/?yp(xo,yo,0) = (ti + xq,— ti + yo,0). 
In particular, for pQ = (xo,yo,0) = (xq,—xq + A;,0) G n the first 
two coordinates of (/>yp(t,Po) are x{t) = t Y xq and y{t) = —t — xq Y k = 
— {t Y Xq) Y k = —x{t) Y k. This concludes the proof. □ 

Proposition 13. The plane ttq is invariant by the flow of Z§, given by (1141) . 

Proof. By Proposition [12] and Proposition [4| we get that ttq is invariant by 
the flow of Zf. □ 
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Remark 5. Since, by Proposition^ we get that all trajectories of Zq con¬ 
verge to TTo we obtain that, for e sufficiently small, the same holds for the 
trajectories of Z§. 

Proposition 14. Consider an integer k > 0 and Z§ given by (1141) . where 
either p = f or p = i. Then z[ has exactly k limit cycles and Z\ has infinite 
many limit cycles, all of then situated in ttq. 

Proof. According to Remark [5l vro is a global attractor for Zf. Also, accord¬ 
ing to Proposition 1131 ttq is invariant by the flow of Zf. So, if there exists 
limit cycles, then they are situated at ttq. Moreover when we restrict the 
flow of Z§ to TTo, by Propositions 1121 and 141 the hxed points of the hrst return 
map ifzp = PYf ° Tx occurs when t = t^ = 4xo. So, take po = (xg, —xg, 0) 
and we get 

(16) TzP(Po) = f>Yf(2xo,-po) = (xo,-Xg,Ff(xg,-Xo)). 

When p = /, by Item (ii) of Lemma [TOl 
ip^p(xg, —xo,0) = (xg, —xg, 0) xg = js with £ > 0 and j = 1, 2,... , A;. 

Therefore, z/ has k limit cycles, all of then situated in ttq. When p = i, by 
Lemma [TT] 

Tzp(xo,-xg,0) = (xo,-xg,0) Xg = jj with j = 1,2,... . 

Therefore, Z\ has infinite many limit cycles, all of then situated in Tig. □ 

Proposition 15. All limit cycles in Proposition are hyperbolic. More¬ 
over, for e > 0, if j = even then it is attractor and if j = odd then it is 
repeller. 

Proof. In fact, in order to prove this we must consider the expression of the 
derivatives in Item (hi) of Lemma uni and Lemma m Observe that when 
A: = 0 in Proposition [T4] there is not limit cycles. In this Item (hi) of Lemma 
do] also is true and the origin is an attractor equilibrium of the system. □ 

Remark 6. Observe that when A; = 0 m Proposition\lf \ the system has not 
limit cycles. In this Item (Hi) of Lemma{T^ also is true and the origin is an 
attractor equilibrium of the system. 

4.2. About the convergence of the trajectories of Zf, given by (fT41) . 

Now we will proceed the analysis of the behavior of the trajectories of (I14p 
in a neighborhood of the T-singularity (at the origin). 

As stated in Propositions 1131 and Remark [S] TTg is invariant by the flow of 
(jldh and all trajectories of (I14h in a neighborhood of the T-singularity con¬ 
verge to it. However, after the perturbation imposed to Zg, the asymptotic 
behavior of the trajectories can drastically changes. The biggest change 
occurs in V~ when p = f and e > 0. Following Proposition [T5l we can 
separate the analysis of this situation, essentially, in three cases: 

• When k = 0 we get that Zf does not presents limit cycles. So, the 
trajectories in becomes increasingly distant from the origin until 
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certain moment. However, since the origin is an attractor when we restrict 
the analysis to vro, there exists a moment from which the trajectory starts 
to converge to origin. So, the T-singularity at the origin of (jl4[) is asymp¬ 
totically stable. 

• When A; 7 ^ 0 is even we get that Z§ presents k hyperbolic nested limit 
cycles and r^, the bigger of them, is attractor. So, the trajectories in 
becomes increasingly distant from the origin until certain moment. However, 
since the bigger limit cycle is an attractor when we restrict the analysis to 
ttq, there exists a moment from which the trajectory starts to converge to 
it. The same holds for all attractor limit cycles at the interior of T^. So, 
the T-singularity at the origin of (jl4p is asymptotically stable. Also, each 
one of the repeller hyperbolic limit cycles at the interior of T^ is repeller 
and each one of the attractor hyperbolic limit cycles is an attractor for the 
trajectories in V~ U TTg. In fact, there are topological half-cylinder of orbits 
converging to each it one of them. 

• When fc 7 ^ 0 is odd we get that Zf presents k hyperbolic nested limit 
cycles and T^, the bigger of them, is repeller. So, there are trajectories in 

leaving any neighborhood of the origin. Also, there are trajectories in 
converging to the origin or, when A; > 3, to an inward hyperbolic limit 
cycle. Moreover, the T-singularity at the origin of (|14p is asymptotically 
stable. Each one of the repeller hyperbolic limit cycles is a repeller and each 
one of the attractor hyperbolic limit cycles is an attractor for the trajectories 
in V~ U ttq. In fact, there are topological half-cylinder of orbits converging 
to each it one of them. 


5. Proof of main results 
Now we prove the main results of the paper: 

Proof of Theorem A:. Item (a): It follows from Remark 01 
Item (b): It follows from Propositions 1141 and 1151 
Item (c): It follows from Propositions [9] and Remark [5l 
Item (d): It follows from Remark [H 

□ 

Before to prove Theorem B, let us define the classical notion of codimen¬ 
sion of vector fields. 

Definition 16. Consider 0(VP) C 14^ a small neighborhood a vector field 
W. We say that W has codimension k if it appears k distinet topological 
types of vector fields in Q{W). 

Proof of Theorem B:. Suppose that the codimension of the T-singularity 
in[2]is m < oo. Then, in a neighborhood of Zq there are PSVFs of, at most, 
m distinct topological types. This is a contradiction due Theorem A. So, 
the codimension of this singularity is infinite. □ 
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6. Conclusion 

As usual, the main aim of the perturbation theory is to approximate a 
given dynamical system by a more familiar one, regarding the former as a 
perturbation of the latter. The problem is to deduce dynamical properties of 
the unperturbed from the perturbed case. In this sense, we focus on certain 
PSVFs Zg which are deformations of Zq. 

In what follows we present a rough description of some properties that 
Zq = {X,Yq) and Zf = (A, T/) enjoy simultaneously. Some of them are: 

• The origin is an equilibrium point of both Zq = (A, Tq) and Zf = 
(A, y/). Moreover, the tangency sets Sx and Syq (resp. Syp), meet 
transversally at the origin, 

• The plane ttq = {y + x = 0} is Zq and Zf-invariant. 

• The sliding vector helds associated to Zq and Zf are topologically 
equivalent. 

Also, there are properties that Zq and Z§ do not enjoy simultaneously. 
For example: 

• The T-singularity at the origin of Z§ is asymptotically stable. 

• The PSVF Zq restricted to tto has a center while in Zf this center is 
perturbed giving rise to hyperbolic limit cycles at ttq. 

Finally, as a conclusion of this paper we are able to say that the T- 
singularity of the PSVF ([2]) has infinite codimension. 

7. Appendix 

Now we illustrate the theoretical analysis performed by means of some nu¬ 
merical simulations. In the next illustrations we use the computer program 
entitled ” Mathematica”. 

We use the following line of commands: 

solutions = Table[First[NDSolve[{x'[t] == > 0, 

-1 - (x[t] + y[t]),l],y'[t] == If[z[t] > 0,1 - {x[t] -hy[t]),-l], 
z'[t] == If[z[t] > 0,-y[t],x[t]],x[0] == 6l,y[0] == -6», 
z[0] == 0}, {x, y, zj, {t, 0, 2.5}]], {0,0.1, 27r - 0.1,0.1}]; 

solutions2 = Table[First[NDSolve[{x'[t] == If[z[t] > 0, 

-1 - {x[t] +y[t]),l],y'[t] == If[z[t] > 0,1 - {x[t] -hy[t]),-l], 

z'[t] == If[z[t] > 0,-y[t],x[t]],x[0] == cos[6l]/4,y[0] == sin[6»]/4, 
z[0] == 0}, {x, y, 4, {t, 0,2.5}]], {0,0.1, 27r - 0.1,0.1}]; 

cl = ContourPlot3D[z, {x, —.4, .4}, {y, —.4, .4}, {z, —.5, .5}, 
Contours —)• 0,Mesh —>• False]] 


c2 = ContourPlot3D[y, {x, —.5, .5}, {y, — A, .4}, {z, —.5, .5} 
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Contours —)• 0,Mesh —)• False]] 


c3 = ContourPlot3D[x + y, {x, —.5, .5}, {y, — A, .4}, {z, —.5, .5}, 
Contours —)• 0,Mesh —)• False]] 

The command 

p3 = ParametricPlot3D[Evaluate[{x[t],y[t], z[t]}/.solutions], {t, 0, 2.5}, 

PlotStyle —)• {Thickness[.0015], Red}]] Show\p3, ImageSize Large] 

generates Figure [H 



The command 

p4 = ParametricPlot3D[Evaluate[{x[t],y[t], z[t]}/.solution.s2], {t, 0, 2.5}, 

PlotStyle — {Thickness[.0015], Red}]] Show[p4:, ImageSize Large] 
generates Figure [71 
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Figure 7. Phase Portrait of trajectories in a neighborhood of the origin. 
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